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In Jaouadi et al. [Phys. Rev. A 83, 023616 (2011)] the authors derive an analytical finite-size 
expansion for the Bose-Einstein condensation critical temperature of an ideal Bose gas in a generic 
power-law trap. In the case of a harmonic trap, this expansion adds higher order terms to the well- 
known first order correction. We point out a delicate point in connection to these results, showing 
that the claims of Jaouadi et al. should be treated with caution. In particular, for a harmonic 
trap, the given expansion yields results that, depending on what is considered to be the critical 
temperature of the finite system, do not generally improve on the established first order correction. 

For some non-harmonic traps, the results differ at first order from other results in the literature. 


The critical temperature for Bose-Einstein condensa¬ 
tion (BEC) of an ideal gas trapped in a generic carte¬ 
sian power-law potential in the thermodynamic limit was 
given in [l|. In the case of power-law traps with spher¬ 
ical or axial symmetry, this temperature was given in 
[1,Q. We will denote it by T° (the superscript 0 meaning 
thermodynamic limit). Following the first realizations of 
BEC with dilute trapped gases, corrections to due to 
the finite number of trapped particles were given, in the 
case of an isotropic harmonic trap, in [^. Shortly after, 
this was generalized to the anisotropic harmonic trap in 
the result reading 


This is eq. (30) in @ (apart from a small difference in 
notation for to a and ujg). We see that this results ex¬ 
tends the older one to higher order. The authors note 
that the second term in ([2]) is of order instead of 

the more expectable which makes the new correc¬ 

tion important for N < 10®. They also present specific 
results for some non-harmonic power-law traps. In this 
comment, we point out that these results depend on how 
the critical temperature is defined for the finite system. 
Specifically, other possible approaches can lead to con¬ 
clusions that differ significantly from the conclusions of 
Q, for some choices of the parameters. 
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where C is the Riemann zeta function, Wa and uig are the 
arithmetic and geometric means of the trap frequencies 
respectively and N is the number of particles. 

Actually, it is well known that for such a system (with 
a finite number of particles) BEC does not strictly exist 
as a sharp, mathematically defined phase transition. In¬ 
stead, it gradually happens over a narrow temperature 
window as the temperature is lowered. It follows that 
eq. © cannot be taken too seriously. It merely provides 
a reference value for the location of this temperature win¬ 
dow (and it is very good at doing that). 

In 0 the authors report an analytical finite-size ex¬ 
pansion for ATc/T®, which is given formally to all or¬ 
ders, in powers of xq = Eo/^ksTc), where Eq is the 
single particle ground level energy. Their framework 
is that of a general anisotropic power-law trap. In 
the specific case of a harmonic trap, we have Xq = 
{3uJa)/{2ujg){(4{‘4)/NY/^{T^/Tc) and truncating at the 
second correction term, their result reads 
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In order to obtain numerical values for the BEC crit- 

(1) ical temperature of a finite system, some specific defini¬ 
tion for this temperature, which works for such systems, 
is necessarily involved. Two common criteria in the liter¬ 
ature are the maximum of the specific heat curve against 
temperature (vanishing first derivative of C{T)) and the 
inflection point of the condensate fraction curve against 
temperature (vanishing second derivative of Ngr/N{T)). 
We will denote the critical temperatures obtained using 
these two criteria by T^ax and Ti„fl, respectively. A mul¬ 
titude of possible criteria and their appropriateness, in¬ 
cluding the two above, was discussed early on in Q, the 
authors concluding in favor of the specific heat maximum 
as the best one. After the first successful BEC experi¬ 
ments with trapped gases, one of those authors @ consid¬ 
ered again the criteria for BEC, this time in the specific 
context of harmonic traps, concluding in favor of eq. o 
as is (without further correction terms), as providing the 
best definition of a critical temperature. Nevertheless, 
we stress that this is a subjective matter. 

The criterium used in @ is implicit in the adopted 
procedure. The authors work in the context of the local 
density approximation. In this setting, the thermody- 
jY-i /2 namic sums over states are approximated by integrals 
in phase space, which are then integrated over the mo- 

(2) mentum, yielding thermodynamic quantities in the form 
of spatial integrals. The critical temperature was then 
obtained by setting the critical value of the chemical po¬ 
tential at fXc = Eq and omitting from the spatial integral 
for N the trap region where the potential is less than 
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FIG. 1. Rescaled temperature Tmax/Tc (solid line), Tinfl/T® 
(long-dashed line), T^/T^ according to eq. ([T} (short-dashed 
line), TcjT^ according to (dotted line), Tx%/T^ (long- 
dashed/dotted line) andTo,i%/T,? (short-dashed/dotted line). 
See text for definitions of these temperatures. The inset 
zooms in on the large N region. It shares a common hori¬ 
zontal axis with the larger plot. 


Eq (the thermodynamic limit critical temperature is re¬ 
covered by setting Eq = Q). The physical significance of 
this criterium for the critical temperature is unclear to 
us (as is does not seem to be based on a specific physi¬ 
cal attribute of the system) and motivated the numerical 
investigation of different possible criteria which we now 
undertake. 

We have computed numerically, directly from the ther¬ 
modynamic sums, Tjjiax and Ting. We also computed two 
additional temperatures, which might serve as reference 
points: the temperatures at which the condensate frac¬ 
tion is 1 % and 0.1%, which we denote by Ti % and 10.1% 
respectively. In Fig. [1] we present the ratio of these four 
temperatures to as functions of the number of par¬ 
ticles for an isotropic harmonic trap. Also shown are 
the usual first order correction, given by eq. o, and the 
higher order correction of given to order in 

eq. ©• For the numerics, we have actually used also 
the three next order terms of this expansion (not shown 
in eq. ([2|)) to assure its convergence. We see that the 
two criteria, specific heat maximum and inflection point 
of condensate fraction, give similar results, with Tc/T^ 
lowering as the system becomes smaller. As for the tem¬ 
peratures at fixed condensate fraction, Ti% and Tq i%, 
as the system becomes smaller they inflect upwards at 
some point. This is a reflection of the fact that as the 
system becomes smaller the phase transition becomes in¬ 
creasingly diluted and the condensate fraction eventually 
becomes significant even above T°. The corrections to 
given by eq. H]) and by the result of follow interme¬ 
diate paths. We see that the ambiguity in what can be 
considered as the transition region is very marked to¬ 
wards the low N range. In particular, except perhaps for 
large N, it does not seem very meaningful to talk about 
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FIG. 2. The same rescaled temperatures as in Fig.[T]are shown 
as functions of the anisotropy parameter s for the case of A = 
10® particles in a disc shaped (upper plot) and cigar shaped 
(lower plot) harmonic trap. The meaning of the several curve 
patterns is the same as in Fig. [T] 


the critical temperature (unless we previously fix it with 
some - necessarily arbitrary - definition). We must point 
out that the authors in Q did not consider values of N 
lower than 10^ (their Fig. 1 has the N range 10^-10®). 

We now consider anisotropic traps. In Fig. [2] we plot 
the same quantities for axially symmetric disc shaped 
and cigar shaped harmonic traps, with anisotropy factor 
s = LOijujj (where i = axial, j = radial for the disc trap 
and the opposite for the cigar trap), as functions of s for 
N = 10®. We see that as the anisotropy increases the 
transition happens at increasingly lower T/T^. However, 
unlike when N is lowered, when s is increased the Ti% 
and Tq 1 % curves do not deviate significantly from the 
Tlnax and Tinfl curves. This is because in this case T^/T^ 
is reduced without diluting the phase transition, which 
remains sharp. We note that for the highly anisotropic 
cigar shaped trap, the critical temperature given by the 
result of @ deviates clearly from the critical temper¬ 
atures given by the other criteria we have considered, 
whereas the one given by eq. o does not. (If the expan¬ 
sion of 0 were truncated at order i.e., eq. the 

respective curve would be slightly above its current po¬ 
sition in both plots of Fig. [21 therefore deviating slightly 
more from the other curves.) 

Summarizing, in the case of a harmonic trap, the re¬ 
sult for ATc/T° given by [1] does not seem to provide 
a better pointer to the transition region than the first 
order result of eq. O- However, due to the lack of a 
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single, well-defined critical temperature, this observation 
is somewhat subjective (except in the case of the highly 
anisotropic disc shaped trap, where the result of [B] falls 
above the transition region). Overall, we believe it is not 
helpful to talk about higher order finite-size correction 
terms to the critical temperature, unless we previously 
fix it by adopting some physical criterium for its defi¬ 
nition. Likewise, the claim in @ that eq. © for the 
critical temperature is innacurate for N < 10®, the new 
term in eq. @ becoming significant in this range, is not 
meaningful on its own, without the respective criterium. 

In what regards non-harmonic power-law traps, we 
limit ourselves to pointing out the existence of other re¬ 
sults in the literature which are in disagreement with the 
result of @ for power-law potentials between cubic and 
power 6 already at first order. Specifically, in Q the first 
order result ATc/T° = —xoC{n)/Kv + 1 )C(^ + 1 )]: where 
we are using the notation of [g, was derived for general 
power-law potentials satisfying ry > 1 . 77 is a constant 
defined in ^ that is related only to the powers of the 
power-law trap. 77 > 1 roughly means a confining poten¬ 
tial less steep than a power-law of power 6 . Eq. © is 
the 77 = 2 particular case of the above formula. The au¬ 
thors in used a high temperature finite-size expansion 
for N which takes into account the discreteness of the 
energy levels. The simpler method of m yields exactly 
the same result. Although in the end result was pre¬ 
sented for the harmonic potential only, the generalization 
to power-law potentials with 77 > 1 is trivial from the in¬ 
termediate results therein. The method of [l0| was used 


11 | in the specific case of a quartic power-law trap 


(for which 77 = 5/4), yielding the 77 = 5/4 instance of the 
same formula. Now, this formula clashes with the result 
of in the cases where 1 < 77 < 3/2 (potential between 
cubic and power 6 ). Indeed, while the first order shift 
ATc/T^ it predicts is linear in xq and negative, the one 
predicted by Q for 77 < 3/2 is of order Xq~^^^, therefore 
more pronounced, and positive. Agreement at first order 
happens only for 77 > 3/2 (roughly, power-laws less steep 
than cubic, including the harmonic trap). 


The disagreement for 1 < 77 < 3/2 might be un¬ 
derstandable on the grounds that the critical tempera¬ 
ture is definition dependent (or ill-defined) for finite sys¬ 
tems. However, note that this means that the conclusions 
drawn in Q for non-harmonic traps also depend on the 
critical temperature definition. In particular, from their 
expansion the authors conclude that in the 77 = 3/2 case 
(e.g., a cubic power-law such as V{r) = ar^) the finite- 
size effects on Tc completely cancel and that for very 
steep power-laws the critical temperature shift can be 
very significant and positive, with possible implications 
for the experimental realization of BEC in such traps. 
These conclusions are contrary to the finite-size correc¬ 
tions obtained in iEl. Thus, they cannot be taken 
in a strict sense and should be considered with caution. 
Moreover, the method of seems to us more adequate, 
in principle, for finite-size calculations than the method 
of [y], since it takes into account the discreteness of the 
energy levels. We believe this particular matter deserves 
further investigation. 
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